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LONG-TERM  GOALS 

The  long-term  goals  of  this  work  are:  (1)  to  develop  a  reverberation  model  for  predicting  the 
reverberation  level  and  the  echo-to-reverberation  ratio  in  shallow  water  (SW),  and  (2)  to  characterize 
the  seabed  scattering  in  a  frequency  band  of  100-3000  Hz  from  long-range  broadband  reverberation 
measurements. 

OBJECTIVES 

The  scientific  objectives  of  this  work  include:  (1)  to  integrate  the  energy  flux  method  for  SW 
reverberation  with  physics-based  seabed  scattering  models  for  reverberation  prediction  and  inversion, 
and  (2)  to  analyze  the  mechanisms  of  low-frequency  (LF)  seabed  scattering,  obtained  from 
reverberation  data  and  the  Biot  seabed  geo-acoustic  model. 

APPROACH 

(1)  The  energy  flux  (angular  power  spectrum)  method  for  reverberation,  based  on  the  W.K.B 
approximation  to  the  normal-mode  solution  [1],  is  integrated  with  physics-based  seabed  scattering 
models  [2]  to  get  SW  reverberation  expressions;  (2)  The  effective  Biot  acoustic  model  for  sandy  and 
sand-silt  mixture  bottoms,  derived  from  broadband  LF  field  measurements  at  many  locations  [3],  is 
used  to  model  two-way  sound  propagation  in  the  reverberation;  (3)  The  broadband  reverberation  data 
are  used  to  characterize  the  seabed  scattering;  (4)  The  LF  seabed  scattering  cross-sections,  derived 
from  the  long-range  reverberation  data,  will  be  used  to  estimate  the  geo-physical  parameters  of  the 
seabed  roughness  and  the  sediment  inhomogeneity,  and  to  analyze  LF  seabed  scattering  mechanisms. 

RESULTS 

(1)  The  energy  flux  method  for  SW  reverberation  is  integrated  with  the  rough  bottom  scattering  (RBS) 
model  and  the  sediment  volume  scattering  (SVS)  model.  This  integration  directly  and  intuitively 
results  in  general  expressions  for  SW  reverberation  in  the  angular  domain  and  in  the  modal  domain. 
The  resultant  reverberation  expressions  can  conveniently  be  used  for  predicting  SW  reverberation  for 
the  Pekeris  waveguide  or  non-Pekeris  waveguide,  respectively.  The  reverberation  intensity  in  shallow 
water,  derived  from  the  energy  flux  method,  can  be  expressed  by  a  summation  over  normal  modes: 
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where  kn  is  the  longitudinal  wave  number  of  the  nth  mode,  and  (3n  is  the  modal  attenuation  coefficient. 
r  is  distance,  and  h  is  the  water  depth;  z0  is  the  source  depth  and  z  is  the  receiver  depth.  aw  is  the 
water  absorption  coefficient.  <j\0m,0n\  is  the  bottom  scattering  cross  section  per  unit  area;  6m(h)  is  the 
incident  grazing  angle  of  mode  m  at  the  bottom;  and  0n  ( h )  is  the  scattering  angle  of  mode  n. 
kw(z )  =  2n f  /  cw(z )  is  the  wave  number  in  the  water.  r0  is  the  signal  pulse  duration.  D(z)  is  a  minor 
modification  factor  for  overcoming  a  problem  when  a  receiver  is  near  a  turning  point, 
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(2) 


For  arbitrary  velocity  profiles,  Eq.  (1)  can  be  used  to  calculate  SW  long-range  reverberation  intensity 
by  using  the  classic  cross  sections  of  seabottom  scattering  and  any  available  normal-mode  code.  This 
only  requires  two  parameters:  real  and  imaginary  parts  of  the  modal  eigenvalue  km  ,  and  /3m  . 


(2)  The  integration  of  the  energy  flux  method  with  seabed  scattering  models  also  results  in  a  simple 
relationship  between  the  classic  boundary  scattering  cross  sections  and  the  modal  scattering  matrix  in 
SW  waveguides: 


MSM  ( m ,  n) 
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Here  MSM(m,ri)  is  defined  as  the  modal  scattering  matrix  in  the  SW  waveguide  between  an  incident 
mode  m  and  a  scattering  mode  n  .  Om(/z)  is  the  amplitude  of  mh  mode  at  the  water-bottom  interface 
where  the  modal  angle  (9m  =  cos  /  kw(h)\  ■  cr(0m,0n)  can  be  either  the  RBS  cross  section  or  the 
SVS  cross  section.  Vww  is  the  sound  reflection  coefficient  with  the  scripts  ww  added  to  indicate  that 
the  incident  and  reflected  fields  are  both  in  the  water. 


(3)  The  data-model  comparisons  show  that  the  reverberation  can  almost  be  equally  well  predicted  by 
using  the  RBS  model  and  SVS  model.  Thus,  the  scattering  parameters  for  the  RBS  model  and  the  SVS 
model  have  been  inverted  from  reverberation  data.  However,  the  reverberation  level  as  a  function  of 
time,RL(t)  ,  at  one  single  frequency  cannot  uniquely  detennine  the  seabed  scattering  parameters.  As 
an  example,  Figure  1  shows  the  internal  parameter  coupling  for  the  RBS  model  between  the  spectral 
exponent  y1  and  the  spectral  strength  a>2 :  three  different  pairs  of  y2  and  <x>2  would  result  in  an 

identical  RL(t) .  An  increase  (decrease)  in  y1  can  be  compensated  for  by  an  increase  (decrease)  in®2 . 


(4)  The  wideband  reverberation  data  as  a  function  of  frequency,  RL(f ) ,  at  given  reverberation  times 
can  uniquely  determine  a  set  of  the  bottom  roughness  spectra  or  the  sediment  volume  scattering  cross 


2 


section.  As  an  example,  Figure  2  shows  that  the  wideband  reverberation  data  are  used  to  uniquely 
invert  a  set  of  seabed  scattering  parameters  for  the  RBS  and  SVS  models.  The  bottom  roughness 
spectrum  and  the  sediment  volume  scattering  cross  section,  inverted  from  the  broadband  long-range 
reverberation  data  at  the  ASIAEX  site,  exhibit  different  characteristics  from  those  commonly  used  in 
HF.  For  the  HF,  the  roughness  spectrum  exponent  is  restricted  to  the  range  of  2  <  <  4  with  a  mean 

value  of  3.0;  the  LF  roughness  spectrum  exponent  from  long-range  reverberation  data  has  smaller 
values  with  a  mean  value  around  y2  «  1.3  .  The  HF  sediment  volume  scattering  cross  section  ov  is 
generally  assumed  to  have  linear  frequency  dependence.  However,  Figure  3  shows  that  the  <jv  , 
inverted  from  long-range  reverberation  data,  exhibits  much  stronger  frequency  dependence  in  the  150- 
2500  Hz  range:  av  =  (0.00238  ±  0.00004)/'3  805±0021) ,  here  f  is  frequency  in  a  unit  of  kHz. 

(5)  All  the  LF  seabed  scattering  parameters  for  the  RBS  model  and  the  SVS  model,  inverted  from  the 
long-range  (3  to  12  s)  reverberation  data,  are  listed  in  Table  I.  The  standard  deviations  between  the 
measured  reverberation  data  and  the  theoretical  predictions  for  all  the  frequencies  are  also  listed  in 
Table  I.  These  seabed  scattering  parameters  for  the  RBS  model  and  the  SVS  model  may  be  used  for 
future  analyzing  which  one  (or  both)  is  dominating  LF  scattering  mechanism  at  low  grazing  angles.  In 
Figure  2  and  Table  I,  a  (in  a  unit  of  dB/rri)  is  the  sound  attenuation  in  the  sea  bottoms  that  obeys  the 
Biot  model  [3], 

IMPACT/APPLICATIONS 

The  resultant  reverberation  expressions  with  the  physics-based  seabed  scattering  models  can 
conveniently  be  used  to  predict  the  SW  reverberation  and  to  characterize  the  LF  seabed  scattering. 
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ASIAEX  site  f  =  2500  Hz  SD  =  50  m  RD  =  56.5  m 


Fig.  1.  Internal  parameter  coupling  between  the  spectral  exponent  y2  and  the  spectral  strength  co2 
for  the  RBS  model.  Three  different  pairs  of  y2  and  co2  result  in  a  identical  RL(t ) .  An  increase  in 

y2  can  be  compensated  for  by  an  increase  in  <x>2 . 
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Fig.  2.  Wideband  reverberation  data  at  the  ASIAEX  site  are  used  to  invert  a  set  of  seabed 
scattering  parameters  for  the  RBS  and  SVS  models. 
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Fig.  3  The  sediment  volume  scattering  cross  section,  inverted  from  long-range  reverberation  data  at 
the  ASIAEX  site,  exhibits  strong  frequency  dependence:  av  =  0.00238 /3'805  (/  in  kHz ) 


Table  I.  The  LF  seabed  scattering  parameters  for  the  RBS  model  and  the  SVS  model,  inverted  from 
long-range  (3  to  12  s)  reverberation  data  at  the  ASIAEX  site 


Models 

Rev.  time  and 
rang-to-depth  ratio  \ 

RBS  :  Bottom  roughness  spectrum 

W(AK)  =  co2/(AKy> 

SVS:  Volume  scattering  cross 
section:  (Jv  =  <72a"' 

Yi 
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